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The CLEO-c detector at CESR has begun to collect large data sets of e + e~ — > cc events 
in the center-of-mass energy range \fs = 3 — 5 GeV that will greatly enhance the size of 
the world's data sets of events produced at charm threshold. Preliminary results for the 
pseudoscalar decay constant f D + and absolute D-meson hadronic branching fractions from 
an initial integrated luminosity C = 57 pb -1 collected at y^i = ^(3770) are described. These 
charm measurements are also relevant for B physics. 



1 Leptonic D + Decays 

Precision measurements of leptonic (and semileptonic) decays in the charm sector are vital for 
determining various elements of the CKM matrix that describes the mixing of quark flavors and 
generations induced by the weak interaction. They also can also be used to relate CKM matrix 
elements to mixing measurements in the BB system and to provide important experimental 
checks of Lattice QCD calculations. Measurements of pseudoscalar decay constants like /^+ are 
of particular interest. 

The lowest order expression for the leptonic branching fraction for the reaction D + — > l + v 
is given by^ 

B{D - Iv) = ^M D+ mf(l - ^-)fl + \V cd \\ D+ , (1) 

where Gp is the Fermi coupling constant, Mp is the D + mass, mi is the mass of the final state 
lepton, f D + is the parameter encapsulating the strong physics of the process, Vd is the CKM 
matrix element encapsulating the weak physics and quantifies the amplitude for c and d quark 



mixing, and Td+ is the D + lifetime. A measurement of B{D + — > /i + z^) coupled with PDG values 
for V c d and Tq+ yields a measurement of /d+. CLEO concentrates on the reaction D + — > u + u 
since for the final state charged lepton u + the branching fraction is relatively large and there 
are fewer final state neutrinos than for the r + case. 

We use 57 pb" 1 of data collected with the CLEO-c detector at = \&(3770). The raw 
data set consists of D°D°, D + D~ and continuum events, with possibly a small admixture of 
t + t~ and two-photon events. Events are divided into a "D- tag side" and a "signal side," where 
the tag side is defined by attempting to fully reconstruct a candidate D~ decay into any one 
of the 5 modes: D~ — > K + -k~tt~, K + tt^tt^tt° , Kgir~, Kgir~ ir~ tt + , Kgir^ir . (Charge conjugate 
modes are implied in this analysis.) Charged kaons and pions are identified by using both 
dE/dx information from the central tracking chamber and information from the ring imaging 
Cherenkov (RICH) counter. Neutral pions are reconstructed using photon shower shape and 
location in the electromagnetic calorimeter and Kg particles are reconstructed from a kinematic 
fit of a pair of charged pions to a displaced vertex. Cutting on the "beam constrained" D~ mass 
ttid = \J E^ eam — (J2Pi) 2 , where the sum is over final state particles, after fully reconstructing a 
D~ candidate, is also a powerful selection tool. The distribution in data for as a function of 
D~ tag is shown in Fig. ^where the fit curves are the superposition of Gaussian signal functions 
and polynomial background functions. We find 28 574 ± 207 tag candidates and 8 765 ± 784 
background candidates. 

Signal D + — > u + u candidates are selected after D~ tag candidates have been identified by 
searching for a single, oppositely charged track, presumed to be a muon. The presence of a 
neutrino is inferred by requiring that the measured value of the missing mass squared (MM 2 ) 
variable be near zero, where MM 2 = (Eb eam — E^) 2 — (—po- — P^+) 2 and Pd- is the three 
momentum of the fully reconstructed D~ . We also require that essentially no (< 250 MeV) 
energy in the calorimeter be unmatched to a track. Fig. |2] shows the Monte Carlo distribution 
of the MM 2 variable for D + —* u + u signal candidates as a function of D tag. The resolution 
is essentially independent of tag mode and confirmed by comparing the resolution of the MM 2 
distribution for both data and Monte Carlo decays of the type D~ — > K s ir~, where the same 
requirements are used as for the u + i> search but with the extra requirement of a single detected 
Kg. Finally, we set our signal search window to be ±2<r in the D + — > fi + u MM 2 data distribution 
with a = 0.028 GeV 2 . 

Fig. El shows the data MM 2 distribution for D + — > fi + u candidate events after having been 
tagged with a fully reconstructed D~ . The insert shows 8 candidates within our search window 
centered about MM 2 = 0. The large peak at positive values of MM 2 is due to D + — > K°ir + 
events. Background events in our search window can arise from non D + — > [x^v decay modes, 
misidentified DD events and from continuum events. The likelihood of these are evaluated using 
Monte Carlo simulation and the results for specific D + decay modes are summarized in Tabled 
Backgrounds from misidentified DD decays and e + e~ — > continuum events are estimated to 
be 0.16 ± 0.16 and 0.17 ± 0.17 events, respectively. The total background is 1.07 ± 1.07 events, 
including systematic errors. 

Table 1: Backgrounds from D + decay modes. 



Mode 


# of events 


D + ^ 

D+ K°7T + 
D+ -» T + V 
D + -» TT°p + U 


0.31 ±0.04 
0.06 ± 0.05 
0.36 ± 0.08 
negligible 



The branching ratio B{D + — > p + v) = N s ig/eN ta g, with the number of background subtracted 
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Figure 1: Distribution in data of the beam constrained mass mo as a function of the D tag mode. The fit curves 
are superpositions of Gaussian signal functions and polynomial background functions, (a) D~ — > K + tt~-k + , (b) 
D~ K+n-Ti+n , (c) D~ -> Kgvr", (d) D~ -> i^7r-7r-7r + , (e) D~ -> T^tt"^ . 
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Figure 2: MM 2 distribution for simulated D + — > fi + v events as a function of the D tag mode. 
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Figure 3: MM 2 distribution in data using D tags plus a single, oppositely charged track with no unmatched 
energy in the calorimeter. The insert shows 8 events inside the signal search window denoted by the arrows. 



signal events N s i g = 6.9 ± 2.8, the Monte Carlo determined single muon detection efficiency 
e = 69.9%, and the number of D T tags N tag = 28 574 ± 207. The results are 



B(D + ii + v) = (3.5 ±1.4 ±0.6) x 10~ 4 and f D + = (201 ± 41 ± 17) MeV . (2) 

These preliminary results are the first statistically unambiguous signal for D + —* [i + v decay, 
consistent with previous claims of observation^ with theoretical models^. CLEO-c intends 
to collect a total integrated luminosity C = 3fb _1 at yfs = ^(3770) in the next year. 

2 Hadronic D-meson Branching Fractions and cr(e + e _ — > \I/(3770) — > DD) 

The absolute branching fractions of the Cabibbo allowed decays D° — > K~n + and D + — > 
K~n + n + (along with Df —* c/>7r + ) are used to normalize the measurements of branching fractions 
for nearly all other D-meson decays. They are also used in extraction of CKM matrix elements 
from b and c quark decay measurements. CLEO has made preliminary measurements of 5 
important D-meson branching fractions D° — > K~ir + , D° — > K~tt + it , D — > K~7t + it + it~ , 
D + —* K~-k + it + and D + —* i^^7r + , as well as the 2 production cross sections <r(e + e~ — > D°D°) 
and cr(e + e _ — > D + D~) at y^s = ^(3770). Charged conjugate particles and decay modes are 
implied throughout this analysis. 

The technique for branching fraction determination follows one first used by the Mark III 
collaboration^ and relies on first fully reconstructing one of the D-mesons ( "single tag" ) in an 
event to tag it as either D°D° or D + D~. Fully reconstructing the second .D-meson ("double 
tag") then allows the absolute branching fraction measurement of either the D or the D + 
independent of the integrated luminosity or the number of DD events produced. For instance, 
if is the total number of DD events produced, then the number of single tag events iVj 

observed decaying via mode i with branching fraction Bi and detected with efficiency is 
Ni = 2N D £)Bi€i. The number of double tag events iVy with the D-mesons decaying via modes 
i and j is then iVy = 2N D fjBiBj6ij for i ^ j and Na = 2N D fjBfeu otherwise. Finally, the 
absolute branching fraction Bi for a given mode i is found by taking the ratio of double tag 
events (iVy) to single tag events (Ni). Bi = (Nij/Nj)(ej/eij) for i / j and Bi = 2(Nu/Nj)(ei/eu) 
otherwise. Including charge conjugate particles and decay modes, we use ten single tag modes 
and thirteen double tag modes. 

Direct use of these expressions for Bi for combining errors and measurements is problematic 
since Ni and Nj are correlated, as are measurements of Bi for different tagging modes j. CLEO 
solves this problem by using a x 2 fitting procedure that simultaneously fits branching fractions 
for all D° and D + decays, and for the number of D°D° and D + D~ pairs produced. Statistical 
and systematic errors, backgrounds, efficiencies and crossfeed for different modes are accounted 
for directly in the fit so that experimentally measured quantities can be accounted for in a 
systematic fashion. 

Event selection is done in a manner similar to that in the fjj+ analysis. Key analysis variables 
for the final selection of D candidates are the energy difference AE = E(D) — Eo, where E(D) 
is the total energy of the particles of the D candidate and Eq is the beam energy, and the 
square of the D candidate mass M 2 (D)c 4 = El - p 2 (D)c 2 , where p(D) is the magnitude of the 
D candidate's three-momentum. The single tag yield in data comes from a binned maximum 
likelihood fit to the M(D) distribution with line shape parameters determined from from both 
Monte Carlo and data. The fit is a sum of Gaussian and Crystal Ball " signal shapes plus an 
ARGUS background function^. D and D distributions for M(D) are fit together with the same 
signal parameters but independent backgrounds. For example, Fig. 0] shows the data M(D) 
distribution for D° — > K~tt + on the left and D + — > K~tt + tt + on the right. 

The double tag yield in data is extracted from an unbinned maximum likelihood fit to the 
two-dimensional distribution of M(D)-M(D) derived from double tag candidate events. Again, 
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Figure 4: M(D) distribution in data for single tag candidates for the decay mode D° — > K~tt + (left) and 
D + — > K~n + iv + (right). The data are dots and the lines are the fit showing the Gaussian and Crystal Ball signal 
functions and ARGUS background function. 
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Figure 5: Projection of the double tag D°-D° candidate masses onto the M(D) axis for the decay mode D° — > 
K~tv + (and its charge conjugate mode) is shown in the left hand plot while the right hand plot shows a similar 
projection of D + -D~ candidate masses for the decay mode D + — > K~-k + tt + (and its charge conjugate mode). 



shape parameters are determined from a combination of Monte Carlo and data events. As an 
example, Fig. El shows the projection of the fit onto the M(D) axis for the case when the D- 
mesons are reconstructed as D° —* K~ir + and D° — > K + ir~ (left-hand side). The right-hand 
side shows a similar projection for the case D + — > K~ir + ir + and D~ — > K + it~tc~ . 

The branching fraction fit for the D-meson branching fractions listed above as well as N DO p 
and N £)+£)- , uses event yields and efficiencies for ten single tag and thirteen double tag modes. 
Crossfeed and backgrounds are accounted for in the fit. The x 2 of the fit is 8.9 for 16 degrees 
of freedom (confidence level = 92%) and the fit results are shown in Table |2J The branching 
fractions are consistent with, but somewhat higher than, the current PDG values. In nearly all 
cases, our statistical error for the branching fractions and their ratios is less than the statistical 
error of the individual measurements comprising the PDG average. Cross sections are extracted 
using our current integrated luminosity C = 57 pb -1 . The statistical and systematic errors of 
these preliminary results should drop markedly when CLEO accumulates its anticipated 3fb _1 
of integrated luminosity. 

Table 2: Fit results for branching fractions and DD pair yields, and inferred production cross sections. First error 
is statistical and second is systematic. 



Parameter 


Fitted/Inferred Value 


B{D U -> K-TT+) 

B(D° -> K~TT + TT°) 

B(D° K-it + it-tt + ) 

B(D + -» K-TT+TT+) 

B{D + -» K° s ir + ) 
N D+D - 


(3.92 ± 0.08 ± 0.23)% 
(14.3 ±0.3 ± 1.0)% 
(8.1 ±0.2 ±0.9)% 
(9.8 ±0.4 ±0.8)% 
(1.61 ±0.08 ±0.15)% 
1.98 x 10 5 
1.48 x 10 5 


a(e + e~ -» ^(3770) -» D [> D°) 
a(e + e- -» ^(3770) -» D+D~) 
a(e + e- -» ^(3770) -» DD) 


3.47 ± 0.07 ± 0.15 nb 
2.59 ±0.11 ±0.11 nb 
6.06 ± 0.13 ± 0.22 nb 
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